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ABSTRACT 

We present measurements of the galaxy cluster X-ray Luminosity Function (XLF) 
from the Wide Angle ROSAT Pointed Survey (WARPS) and quantify its evolution. 
WARPS is a serendipitous survey of the central region of ROSAT pointed observations 
and was carried out in two phases (WARPS-I and WARPS-II). The results here are 
based on a final sample of 124 clusters, complete above a flux limit of 6.5 x 10 -14 erg 
cm~ 2 s _1 , with members out to redshift z ~ 1.05, and a sky coverage of 70.9 deg 2 . We 
find significant evidence for negative evolution of the XLF, which complements the 
majority of X-ray cluster surveys. To quantify the suggested evolution, we perform a 
maximum likelihood analysis and conclude that the evolution is driven by a decreasing 
number density of high luminosity clusters with redshift, while the bulk of the cluster 
population remains nearly unchanged out to redshift z ss 1.1, as expected in a low 
density Universe. The results are found to be insensitive to a variety of sources of 
systematic uncertainty that affect the measurement of the XLF and determination 
of the survey selection function. We perform a Bayesian analysis of the XLF to fully 
account for uncertainties in the local XLF on the measured evolution, and find that 
the detected evolution remains significant at the 95% level. We observe a significant 
excess of clusters in the WARPS at 0.1 < z < 0.3 and Lx ~ 2 x 10 43 erg s _1 compared 
with the reference low-redshift XLF, or our Bayesian fit to the WARPS data. We find 
that the excess cannot be explained by sample variance, or Eddington bias, and is 
unlikely to be due to problems with the survey selection function. 

Key words: cosmology: observations - galaxies: clusters: general - X-rays: galaxies: 
clusters 



1 INTRODUCTION 

Evolutionary properties of gravitationally bound objects in 
the universe are described by models of structure formation. 
The currently favoured cosmology (flat ACDM) predicts lit- 
tle change in the abundance of galaxy clusters at late times 
when the energy density of the universe becomes dominated 
by Oa. The evolution of cluster abundance depends on the 
growth rate /, which is mainly sensitive to the mean cosmic 
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matter density f2 m as /(«) ~ f2 m (2) 7 , where 7 ~ 0.6 for 
General Relativity (Lindcr 20051. 



Galaxy clusters, the largest objects to have decoupled 
from the Hubble expansion, are particularly interesting for 
studying these properties as a result of their X-ray bright- 
ness. The X-ray emission is the result of bremsstrahlung 
emitted by the hot intracluster medium (10 7 — 10 8 K) which 
contributes more than 80% of the baryonic content of the 
cluster. Therefore, the mass of a cluster can be obtained from 
its luminosity with the use of scaling relations and some sim- 
plifying assumptions ( Kaiser]T986 ) . The X-ray emitting gas 
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has enabled cluster detections out to high redshift (z > 1). 
Hence, X-ray galaxy cluster surveys potentially cover a sig- 
nificant portion of the evolution history of clusters and have 
high statistical completeness, thus providing the leverage 



formed by Mullis et al. ( 2004 ) using the 160SD catalogue 



to place tight cosmological constraints (e.g. Vikhlinin et al. 
|2l^|Mantz et al.|2010) . 

Early predictions of evolution in the number density of 
clusters, e.g. |Kaiscr ( 1986]) , pointed towards strong positive 
evolution - an increase in the number density of clusters 
with redshift. This prediction assumes a matter power spec- 
trum with a power-law form, and that the heating of gas is 
solely by adiabatic compression during the collapse of dark 
matter halos. The first opportunity to test these predictions 
came with the Einstein Medium Sensitivity Survey (EMSS 
Gioia et al. 19901, which detected clusters out to z ~ 0.8. 



Contrary to the theoretical prediction, the first teams to 
test for evolution in the XLF found strong negative evolu- 



tion Gioia et al.|1990||Edge et al.|1990||Henry et al.|1992 | 

These controversial findings heated the debate and together 
with the launch of the ROSAT X-ray observatory gave rise 
to a flurry of attempts to measure evolution in the XLF, 
with some later analyses raising concerns over the Einstein 
results (e.g. |Ellis fc Jones||2002"l ) . 

ROSAT performed an all-sky survey which was used 
to construct large flux limited cluster samples, from which 
the local cluster XLF was accurately determined. There are 



three such surveys: the Brightest Cluster Sample (Ebeling 
et~aLlp98l [2000), the ROSAT All-Sky Survey 1 Brightest 
Sample ( |de branch et al.|1999[ ), and the ROSAT- ESO Flux- 
Limited X-ray (REFLEX) galaxy cluster survey |Bohringer 
et al.|200lj ). These local XLFs act as the crucial baseline for 
quantifying evolution in deeper surveys. 

Once the ROSAT all-sky survey was completed the ob- 
servatory remained available for pointed observations, which 
has resulted in an extensive archive of deep observations, 
providing the ingredients for many serendipitous X-ray clus- 
ter surveys. This includes the Wide Angle ROSAT Pointed 
Surve y ( jScharf et al.|1997||Perlman et al.|2002l|Horner et al. 
2008), the subject of this paper. Similar surveys that probe 



the X-ray universe out to high redshift include the ROSAT 



International X-ray/Optical Survey (RIXOS Castander et 
al.|1995| |Mason et al.|2000|>, the ROS AT Deep Cluster Sur- 
vey (RDCS |Rosati et al.||19"95| |1998| l, the Bright Serendip- 
itous Hi gh-Redshift Archiva l ROSAT Cluster (BSHARC) 
Survey ( jRomer et al . 20001, the Massive Cluster Survey 



(MACS) libeling et al 



2001 ), the Brera Multi-scale Wavelet 



ROSAT HRI (B MW-HRI) survey ( |Moretti et al.|2001||Pan^ 
era et al.||2003]) ROSAT North Ecli ptic Pole ( NEP) Survey 
pienry et al.|2001||GioTa et al.|2001[ ), SSHA RC ([Burke et al. 



1998 



2003), and the 160 Square Degree (160SD Vikhlinin et al 



^ Mullis et al.||2003[|2004|), extended to the 400 Square 
Degree (400SD) survey ( |Burenin et al.|2007| . 

XMM-Newton archival data is also used for surveys 
based on serendipitous cluster detections. Currently in 



progress are the XMM Cluster Survey (XCS Mehrtens et al. 
20TT|), the XMM-New ton Distant Cluster Project (XDCP 
Fassbender et al.|201lj), and the XMM-N ewton eXtra eXtra 
Large (XXL) Survey ( Pierre et al.|201T l. One serendipitous 
galaxy cluster survey is based on Chandra archival data and 
is part of the Chandra Multiwavelength Project (CHaMP 
Barkhouse et al.|2006 l. 

The most recent determination of the XLF was per- 



This work found significant evidence for negative evolution 
of the XLF at the bright end. That is, the number density 
of high luminosity clusters was lower at 0.6 < z < 0.8 than 
in the local Universe. 

In this paper we investigate the evolution of the XLF 
of a sample of 124 WARPS galaxy clusters detected above a 
flux limit of 6.5 x 10 -14 erg cm -2 s _1 over a total area of 70.9 
deg 2 , and covering a wide redshift range (0.02 < z < 1.10). 



The survey design was outlined in (Scharf et al. 19971 



and the catalogues are presented in two separate papers: 



WARPS-I (Perlman et al. 2002) and WARPS-II (Horner 



et al. 2008). The evolution of the WARPS galaxy clusters 



has previously been investigated using phase-I of the sur- 
vey. Jones et al. ( 1998 1 found no significant evolution in the 



logiV — logS relation from the WARPS-I sample and a pre- 
liminary measurement of the XLF (constructed when the 
survey was complete for z < 0.85) was also found to be 
consistent with no evolution ( | Jones et al.|2000 1 . 

This work represents a useful cross-check and extension 



of the Mullis et al. ( 2004 1 results. While the WARPS survey 



covers a smaller area, it is deeper; the 160SD XLF extends 
to z ~ 0.7. Importantly, while both surveys are drawn from 
ROSAT pointed observations, the cluster detection and con- 
firmation stategies differ significantly, allowing us to assess 
the sensitivity of the evolution results to those factors. 

The current paper is organised as follows: §2 briefly re- 
views the WARPS survey and the combined WARPS-I + 
WARPS-II sample. The selection function of the full survey 
is presented for the first time. In §3 the X-ray Luminos- 
ity function is presented for different redshift ranges. Sub- 
sequently, in §4 a maximum likelihood analysis robustly as- 
sesses evolution in the XLF. In §5 and §6 we discuss our re- 
sults and summarise our conclusions. Throughout the paper 
errors are quoted at the 68% confidence level and a ACDM 
cosmology of Ho = 70/i7o km s _1 Mpc -1 , fi m = 0.3, and 
Qa = 0.7, is adopted. All fluxes are corrected for absorption, 
and are quoted in the observer's frame 0.5 — 2 keV band. Lu- 
minosities are converted to the rest frame 0.5 — 2 keV band 
of each cluster. 



2 THE WARPS CLUSTER SAMPLE 

The survey is based entirely on serendipitous detections in 
ROSAT images from pointed observations with the Position 
Sensitive Proportional Counter (PSPC) instrument. Here we 
summarise the key facts and direct the reader to |Scharf et| 
al. ( 1997 1 for full details of the survey methodology. 



PSPC fields were selected based on the following crite- 
ria. The fields are at a Galactic latitude of |&| ^ 20°, have 
exposure times of t Gxp 8ks, are non-overlapping, and the 
original target is not a galaxy cluster or some other source 
such as a bright star that would hamper optical follow-up. 
Out of the ~7000 fields in the HEASARC archive 381 satisfy 
the criteria. 

Sources were detected with Voronoi Tessellation and 



Percolation (VTP) (Ebeling fc Wiedenmann 19931 in an 



annulus of inner radius 3 arcminutes and outer radius 15 
arcminutes. VTP does not discriminate against shape or 
size and is particularly sensitive to sources of low surface 
brightness. WARPS has assessed the efficacy of VTP as a 
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Figure 1. The WARPS clusters (hollow points) as points in the 
Lx-2 plane. Also plotted is the REFLEX sample (solid points), 
which provides the low-redshift reference XLF for this study 
< |Bohringer et al.|200i) . 



source detection algorithm by optically imaging all X-ray 
candidates in WARPS-I lacking counterparts on existing sky 
survey plates. 

All sources with a detection flux of 3.5 x 10~ 14 erg cm -2 
s _1 are detected at 3a by the VTP algorithm. After ex- 
trapolating the detected flux to infinite radius assuming a 
/? profile, a statistically complete sample was defined, com- 
prising 124 clusters above a uniform flux of 6.5 x 10 -14 erg 
cm~ 2 s _1 . Spectroscopic redshifts were obtained for all clus- 
ters, with 2 or more concordant redshifts required to confirm 
a cluster. WARPS did not obtain near-infrared imaging of 
cluster candidates, placing an upper limit on the redshift 
out to which clusters can be detected. This limit is ~ 1.1, 
and the uncertainty arising from this is addressed in section 

In combining WARPS-I and WARPS-II catalogues, 
it was found that background levels were missing for 1 
WARPS-I field and 27 WARPS-II fields. The background 
level of each field is required in order to compute the se- 
lection function, and so these were remeasured from the 
archived PSPC data. The ROSAT PSPC data have been 
reprocessed since the cluster detection was performed, so 
we checked the background measurements for all WARPS-I 
fields using the currently available PSPC data against our 
original measurements. The new measurements were found 
to be ~ 7% lower on average, depending somewhat on the 
source detection algorithm used to exclude sources in each 
field. We thus renormalised the background measurements 
for the 28 missing fields in the combined WARPS catalogue 
by this factor, for consistency with the data used for cluster 
detection. We investigated the impact of this systematic ef- 
fect on the selection function, and found it to be insensitive 
to whether or not this background scaling was applied to 
the 28 missing fields. This is not surprising given the small 
magnitude of the correction and the small fraction of fields 
affected. 

In figure [T] the WARPS clusters are plotted in the 



Figure 2. Survey coverage for three classes of objects as a func- 
tion of redshift. The objects are defined as: elliptical galaxies with 
Lx(0.5 — 2.0 keV) = 1 X 10 42 erg s — 1 and effective core radius 
r c = 50 kpc, groups with L x (0.5 - 2.0 keV) = 1 X 10 43 erg 
s — 1 and effective core radius r c = 100 kpc, and clusters with 
L x (0.5 - 2.0 keV) = 5 X 10 44 erg s" 1 and effective core radius 
r c = 250 kpc. The grey line represents the approximate upper 
redshift limit imposed by the lack of near infra-red follow-up of 
cluster candidates. 



luminosity-redshift plane. The dashed line corresponds to 
the flux limit of 6.5 x 10" 14 erg cm -2 s~\ The fluxes of the 
clusters have been K-corrected to the cluster rest frame as- 



suming an APEC thermal plasma model ( Smith et al. 2001 1, 
for which we set the metallicity to 0.3 Zq. The plasma tem- 
perature required for this conversion was estimated itera- 
tively from the X-ray luminosity using the luminosity tem- 



perature scaling relation of Markevitch ( 1998), although the 



magnitude of the K-correction was insensitive to this choice. 

The selection function for WARPS-I and WARPS-II 
combined is shown in figure [2] and is based on 381 PSPC 
fields. The effective sky coverage for an object of given lu- 
minosity and extent is determined by the performance of 
VTP, the degrading PSF with off-axis angle, and the back- 
ground levels and exposure times of the fields. The variance 
in the field properties of the survey alters the steepness of 
the decrease; e.g. if all the fields were the same, we expect a 
much more sudden drop from 100% to 0%. We find the curve 
for the full survey to be very similar to that for WARPS-I, 



figure 9 in Scharf et al. (1997), indicating that WARPS-I is 



a reliable subset of the full survey. 



3 THE X-RAY LUMINOSITY FUNCTION 

The X-ray Luminosity Function (XLF), conventionally given 
the symbol (j>, is the comoving number density n of objects 
per luminosity interval: 



4>(L 



x.z 



dn.(Lx,z) 
dLx ' 



(1) 
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The Schechter function (Schechter 19761 is the canonical, 



parametric representation of the luminosity function: 



<j>(Lx, z)dLx = (t>* 



l * 



exp 



Lx 
L 



dLx 
L 



(2) 



J x/ V ^x 

where the parameter <ji* normalises the XLF, and a deter- 
mines the steepness at Lx < L x . 

The conventional method to compute the differential 
XLF is the 1/V max method ( |Schmidt||1968| |Avni fc Bahcall| 
1980 1, where Knax denotes the maximum co-moving volume, 



given by 



Vmax — 



tt(fx,re) v ' dz, 
dz 



(3) 



where Q(fx,rg) is the sky coverage as a function of flux 
fx(Lx, z) and angular extent rg(r c ,z) (here r c is the core 
radius of the cluster surface brightness distribution, conven- 
tionally parameterised with a /3-model), and dV(z)/dz is the 
differential, co-moving volume, which is strongly sensitive to 
the cosmological framework. The maximum co-moving vol- 
ume is calculated for all iV galaxy clusters. The XLF is then 
obtained by summing the corresponding density contribu- 
tions per luminosity bin, that is 



E 



(4) 



where the subscript j denotes the j-ih bin. Due to the sen- 
sitivity to the choice of binning, the method is less ideal for 
quantifying evolution. However, it is a conventional way of 
presenting a sample of objects, allowing for visual compari- 
son with previous surveys. 



Alternatively, |Page fc Carrera ( 2000[ ) provide an esti- 
mate of (f>, which expression is obtained by integrating |l]) 
and noting that <f) changes little compared to the survey vol- 
ume element in the volume - luminosity plane, such that it 
can be taken out of the integral, giving 



<t>{Ly 



r L x,t 

JLyr „ 



N 3 

n(fx,r e ) d ^dzdLx 



(5) 



where Lx is the bin centre and Nj is the number of clusters 
in the j-ih bin. 



We apply the method of Page & Carrera ( 2000 1 to ac- 



count for the flux limit of the survey to effectively decrease 
the width of some of the bins, enhancing the XLF. The Page- 



Carrera estimator was also deployed by Mullis et al. (20041, 



who found a marginal increase at the faint end of the XLF 
compared to the Knax estimator. We can complement these 
findings and, since other uncertainties at the faint end are 
significantly larger - in particular those arising from small 
number statistics and fx and re , we conclude that the choice 
of method is unlikely to affect our results and conclusions. 




O WARPS (0.02 <z< 0.3) 
REFLEX (z < 0.3) 



1 Q 

L y (0.5-2.0 keV) [h^ 10 44 ergs _1 ] 



Figure 3. The XLF from the local WARPS sample along with 
the best fit Schechter function of the REFLEX sample. 



in each luminosity bin are provided by Gehrels ( 1986 1, which 



are much larger than the flux measurement errors. 

Good knowledge of the local XLF is essential for study- 
ing its evolution and is provided with great accuracy by the 
ROSAT all-sky survey. The XLF of the local WARPS sample 
of 67 0.02 < z < 0.3 clusters is shown in figure |3] The lower 
redshift limit is set to z = 0.02 below which many clusters 
become too extended relative to the size of the PSPC fields 
to be detected. Over this redshift range the WARPS XLF 
agrees remarkably well with the all-sky samples, represented 
by the REFLEX model in figure [3] 

There appears to be a high number density of clus- 
ters at ix(o.5-2.okcV) ~ 1-5 x 10 43 erg s _1 compared to the 
Schechter function. We note that this feature is also present 



in the local XLF of the 160SD sample (Mullis et al. 2004 



figure 4) We test the significance of this excess in section 
|3.2| and discuss possible interpretations in section [472] 

In figure [4] we show the intermediate and high redshift 
XLFs along with the local REFLEX Schechter function. The 
majority of data points of both the intermediate and high 
redshift XLF are slightly low compared to the local baseline. 
This is a first indication from the data of negative evolution. 
Whether this is significant will be addressed in the next 
section. 



3.2 Expected Versus Observed Numbers 

The expected number of objects in the luminosity-redshift 
plane is obtained by integrating equation |l]) 



3.1 The WARPS XLF 

In order to present the binned WARPS XLF, we divide 
the clusters according to their redshifts to study the local 
(0.02 < z < 0.3), intermediate redshift (0.3 < z < 0. 6), and 
high redshift (0.6 < z < 1.1) populations, similar to Mullis 



et al. (20041. We apply the same Lx binning as Mullis et al 



1 2004) to allow for comparison. Poisson errors on the counts 



0(L X , z)0(/ X , r e ) dzdix- 
dz 



(6) 

As mentioned in section [3] the XLF changes little compared 
to the volume element. Hence we can predict the number 
of clusters for any of the WARPS subsets based on the lo- 
cal reference XLF <^i oca i, the observed XLF for the subset 
(/"observed, and the number of clusters observed in that sub- 
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Figure 4. The XLF from the intermediate and high redshift 
WARPS samples along with the best fit Schcchter function of 
the REFLEX sample. 



set N ob 



■^V observed ^ 



01ocal 



(7) 



If the local reference XLF is a good description of the 
WARPS XLF, and there is no evolution, then AT ex p should 
be consistent with JVobserved for all subsets. 

Using the REFLEX best fit Schechter function as the 
local reference, we compute the expected cluster numbers 
for each luminosity bin in each of the WARPS subsets. The 
results are plotted in figure [5] When integrated over the full 
range of luminosities, 60 clusters are expected from equa- 
tion Q for the low-z subset, instead of the 67 observed. For 
the intermediate-z subset, 67 are predicted instead of the 44 
observed, and for the high redshift subset, the local relation 
predicts 36 clusters instead of the 13 that are observed. The 
differences for the low and intermediate redshift subsets are 
not strongly significant, but the lack of high-z clusters com- 
pared to the local prediction is significant at > 4a, assuming 
Poisson errors on both numbers. 

Figure [5] also illustrates the excess of clusters around 



J X(0.5-2.0kcV) 



1.5 x 10 43 erg s 1 in the low-z subset seen 



in figure [3] Over the two bins with excess counts, there are 
28 clusters observed, while only 14 are predicted by the RE- 
FLEX XLF. This is a significant excess; the probability of 
observing N > 27 for a Poisson distribution with a mean 
of 14 is 6.4 x 10" 4 . The same analysis was also applied to 
the larger low-z subset of the 160SD sample. According to 
equation Q, the number of clusters predicted by the local 
REFLEX XLF over the two bins with excess counts for the 
160SD sample is 18, significantly lower than the observed 
number of 40 clusters [P(N > 39) = 5.3 x 10~ 6 for a Pois- 
son distribution with mean 18]. 



3.3 Evolving Schechter function 

Here we deploy the maximum likelihood analysis first set out 



by Marshall et al. (19831, which fits an evolving Schechter 
function to the distribution of objects in luminosity redshift 
space. The treatment is free from arbitrary binning and with 






0.3<z<0.6 









L (0.5 - 2.0 keV) [hf 10" erg s~ 1 ] 




L,, (0.5 - 2.0 keV) [h^ 10" erg s~ 1 ] 



Figure 5. Expected cluster numbers (solid line) versus observed 
(dots) per luminosity bin for the local, intermediate, and high 
redshift samples. The expected number of clusters is calculated 
per luminosity bin from equation 161. 



the generalisation of Mullis et al. (20041 accounts for flux 
uncertainties. We briefly summarise the method and apply 
it to the WARPS sample. 

The XLF is characterised as an evolving Schechter func- 
tion 



</>(L x , z)dLx = 4>*{z) 



exp 



Lx 



d£ x 



(8) 

The parameters, except for a, are allowed to evolve as fol- 
lows 



4>*{. z ) = <t>o 



1 + 2 



l + z (Lx) 

1 + z 
1 + Zo(L% 



(9) 
(10) 



where 4>o an d £ x ,o are adopted from the local XLF. Due 
to the flux limit of the surveys, the median redshift zq in- 
creases with luminosity bin and is given by the local XLF. 
A deviation from A — B — indicates evolution. 

To be free from arbitrary binning, the luminosity red- 
shift grid is chosen to be sufficiently fine for there to be either 
1 or clusters in each cell. We achieve this with dz = 0.01 
and dL x = 0.1 x 10 43 /i" 2 erg s^ 1 for the WARPS sample. 
In each cell the expected number of clusters is calculated 

X(L X , z)dLxdz = 4>{Lx, zMfx, r e )^^-dLxdz (11) 

The likelihood function £ describes the joint probability of 
detecting 1 cluster at each occupied cell i and in each 
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Figure 6. Likelihood contours for the evolution parameters A 
and B based on a comparison of the local REFLEX Schechter 
function and the complete WARPS distribution of objects in lu- 
minosity redshift space. Contours show the la, 2a, and 3<r confi- 
dence limits. 




Figure 7. Likelihood contours for the evolution parameters A 
and B based on a comparison of the local REFLEX Schechter 
function and the z > 0.3 WARPS clusters. Contours show the 
la, 2a, and 3a confidence limits. 



empty cell j and is given by 

C = ]^[A(Lx,i, Zi)dLxdz exp [— A(L x ,i, Zi)dLxdz] 

i 

x Y\ CX P [-AfXx.i, Zi)dL x dz] , 



(12) 



which makes use of the Poisson distribution and is valid 
when the number of expected clusters < 1 as expected for 
small cells. The objects are smoothed by a Gaussian in the 
luminosity direction according to the la flux errors. Redshift 
errors are not taken into account, since they are typically 
much smaller. 

We calculate AS = S(A, B) - S{A hast ,B hcst ), where 
S = -21n£. In figure [6] we plot contours of AS = 2.30, 6.17, 
and 11.8, which correspond to the la, 2a, and 3<r confidence 
limits. 

The contours for WARPS are shown in figure [6] with 
evolution measured relative to the local XLF from REFLEX, 
for which we used all 124 WARPS clusters. We find evi- 
dence for negative evolution that is significant at 2a, with 
A = -1.88 ± 0.62 and B = -1.76 ± 0.53. Stronger evidence 
for evolution was measured using only the WARPS clusters 
at z > 0.3. The same methodology was applied, and the re- 
sulting confidence contours on A and B are shown in figure 
[7] The best fitting parameters for both samples are given in 
table [2 



4 DISCUSSION 

4.1 Evolution in the XLF 

The comparison of expected and observed cluster number 
counts, and the maximum likelihood analysis of the un- 
binned cluster population both strongly support negative 
evolution of the XLF. The evolution in (f>* and L* x is degen- 
erate, as is apparent in figure [5J but the net effect is signifi- 
cant, and consistent with decrease in the number density of 




1e+43 1e+44 

L x (0.5-2.0 keV) Id; 10"erg s" 1 



Figure 8. Schechter functions with best-fitting evolution param- 
eters from the maximum likelihood analysis to the full sample. 
The lines show the form of the Schechter function at redshifts 0, 
0.5 and 1. 



massive, high luminosity clusters with redshift, as expected 
in a ACDM hierarchical Universe. This is illustrated in fig- 
ure [8] which shows the best fitting evolution of the REFLEX 
z = Schechter function. 

In the maximum likelihood analysis of section [3~3| we as- 
sume no evolution in the parameter a. Here we test whether 
this is justified by calculating AS for the parameter C de- 
fined as 



l + z 



1 + z (L K ) 



(13) 



and use our best fit parameters A and B from section |373| We 
find Cbest = 0.05 ±0.17, and for A = B = we have Cbost = 
—0.05 ± 0.17. This is reassuring, since in the hierarchical 
picture of structure formation we expect evolution to occur 
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Dataset 


rcdshift 




a 


L x 


A 


B 






10- 7 ft? Mpc- 3 




10 44 ftf 2 






REFLEX 


z < 0.3 


2.94 ± 0.82 


1.690 ±0.045 


2.64 ±0.29 






WARPS ML 


2 > 0.02 


2.94 


1.690 


2.64 


1.88 ±0.62 


-1.76 ±0.53 


WARPS ML 


z > 0.3 


2.94 


1.690 


2.64 


3.60 ± 0.95 


-3.37 ± 0.56 


WARPS Bayesian 


z > 0.02 


3.68 ±0.87 


1.79 ±0.04 


2.59 ± 0.35 


-0.09 ± 1.19 


-0.93 ±0.58 



Table 1. Best fitting XLF parameters. The REFLEX parameters are taken from (Bohringcr ct al. 2002), for a ACDM cosmology. The 
maximum likelihood (ML) fits assumed XLF shape parameters fixed at the REFLEX best fit values. The Bayesian fit used the REFLEX 
values as priors, as discussed in section [4.1. 3| 



at the bright end, whereas a determines the slope at fainter 
luminosities. 

The negative evolution measured in the WARPS XLF 
is in very good agreement with that found for the 160SD 
(particularly for the z > 0.3 subset), and agrees quali- 
tatively with previous measurements of negative evolution 
in the XLF (see figure 8 in |Mullis et al.||2004[). We note 



that the 160SD evolution was measured with respect to the 
BCS XLF, whereas we used the REFLEX XLF for our low- 
redshift baseline, so the agreement in evolution measures 
suggests that the choice of baseline does not strongly affect 
the measured evolution. The uncertainty on the local XLF 
and other systematics affecting the evolution measurements 
are discussed below. 



4-1.1 Uncertainty on z max 

The upper redshift limit z max for the high redshift sample is 
uncertain due to lack of near infra-red follow-up of cluster 
candidates. To our knowledge it is the first time an issue for 
determining the cluster XLF. The RDCS survey included 
near infra-red imaging which resulted in the successful de- 
tection of 4 clusters beyond redshift 1. However, these fall 
below the flux limit of their statistically complete sample, 
which has z ma x = 0.83. The detection and optical confirma- 
tion in WARPS of a cluster at z = 1.028 is consistent with 
our estimate of z max ~ 1.1, but the exact limit will depend 
on the characteristics of the galaxy populations and the pho- 
tometric limits of the optical imaging for each cluster, so is 
not well defined. 

As can be seen in figure [2j for very luminous clusters 
the sky coverage of 100% is maintained well beyond the limit 
set by the optical observations. This means that the prod- 
uct f2(/x,re) x is nonzero in equation (|5j> . Hence, an 
increase in z max suppresses the XLF in those bins that rep- 
resent sufficiently high X-ray luminosities. Although the in- 
fluence of the choice of zi max on the XLF is suppressed by 
the flux limit, a too high value for z max could falsely sug- 
gest negative evolution, whereas a too low value boosts the 
bright end towards the opposite conclusion. 

We tested the robustness of the evolution measurement 
to z max by reducing z max to z = 1 and excluding the 
z = 1.028 cluster from the analysis. The contours in the 
A — B plane were changed negligibly. This is a conservative 
approach, and shows our measured evolution is insensitive 
to the choice of z max . If the effective Zmax of the survey 
is actually larger than the assumed z — 1.1, then our non- 
detection of clusters beyond z = 1.03 would imply stronger 
evolution than that measured here. 



4-1.2 Cluster Surface Brightness 

Although we refer to the statistically complete WARPS sam- 
ple as "flux-limited", in practice it is the X-ray surface 
brightness and not the flux that determines whether or not 
a cluster is detected. To a first approximation, the surface 
brightness is related to the cluster flux by the core radius r c , 
that sets the spatial scale of the surface brightness distribu- 
tion. This then enters the XLF through the computation of 
the detection volumes of the clusters. Given the relatively 
large PSF of the PSPC at the off axis angles considered in 
WARPS, we do not expect our results to be sensitive to the 
choice of core radius, with the strongest effects expected at 
the faint end of the local XLF, which is most sensitive to 
uncertainties in the selection function. 

For each WARPS cluster, a core radius was estimated 
from the PSPC data as the radius at which the surface 
brightness, fitted by a /3-model with ft = 2/3, is a factor 
2 1///3 lower than the central value. The uncertainties on the 
individual core radii are large, but the average r c for WARPS 
is ~ lOOkpc, whereas Chandra observations of clusters show 
an average core radius of ~ 150 kpc (Mau ghan et al.[2 008). 

A disadvantage of the Page-Carrera technique is that 
information about the core radius of the individual clusters 
is difficult to include. Thus when applying this technique, a 
fixed core radius of 102 kpc (average WARPS) was assumed 
for each cluster. However, using the Kn ax technique it was 
possible to investigate the effect of varying r c . We found 
that the Vm ax technique yields nearly identical XLFs for a 
uniform core radius of 100 kpc (average WARPS), 150 kpc 
(average Chandra), and the individual core radii measured 
from the PSPC data. 



4-1.3 Uncertainties on the Local XLF 

In order to assess the impact of uncertainties on the form 
of the low redshift XLF on the measured evolution, we 
adopted a Bayesian approach to fitting the XLF. The poste- 
rior probability distribution for the set of model parameters 
8 = (<j>* , Lx, ct, A, B) given the observed data D (the lumi- 
nosity and redshift of each cluster) is given in the normal 
way by 



P(0|D) oc P(D|0)P(0). 



(14) 



Here the first term on the right is the likelihood function, 
and the second term is the prior probability distribution of 
the model parameters. This approach allows us to adopt the 
REFLEX low-z XLF parameters and their uncertainties as 
priors on <j)*,L x ,a, which can then be marginalised over. 
We adopt weak priors on A and B, simply assigning each a 



© 2012 RAS, MNRAS 000, 000-000 



8 L. A. Koens et. al. 



Gaussian distribution with mean zero, and variance 16. We 



also generalise the likelihood expression from equation (121 
to include the statistical scatter of the measured luminosi- 
ties. This accounts for the possibility that clusters that are 
nominally below the flux limit may be observed to be above 
the flux limit due to our noisy measurement of Lx (this is 
a source of Eddington bias and is discussed further below). 

We divide the Lx, z parameter space into cells i,j with 
coordinates {Lx,i, Zj) and widths (dLx,i, dzj). As before, the 
XLF model predicts a number of clusters in cell i,j as 

Nm 0d ,ij = \(L x ,i,Zj,0)dLx,idzj. (15) 

However, the final number of clusters expected in cell i,j 
includes contributions from all of the other cells in the Lx 
direction, due to the noisy measurement of Lx- The contri- 
bution from a cell at Lx,k, Zj to the number counts in a cell 
at L x ,i,Zj is 



N mo a,ikP(Lx,j\Lx,k,o'k)dLx, : 



(16) 



^exp,ij k 

The probability term here models the measurement noise 
on a cluster with "true" luminosity Lx.k as a Gaussian with 
a mean Lx,k and standard deviation Ok- We model the in- 
creasing precision of the luminosity measurement with clus- 
ter flux by setting Ok to be inversely proportional to the 
square root of the flux at Lx,k,Zj, as expected for mea- 
surements dominated by counting statistics. The constant 
of proportionality is set to give a 15% luminosity error at 
the flux limit, in agreement with the observed clusters. 

The final expected number of clusters in cell Lx.i, Zj is 

then 



/] A / ' mo d,ifc-P(ix,i|-Lx,fc, Vk)dLx,j 



and the likelihood function is then 

P(P\0) = ]JP(N ohs , lJ \N c: 



(17) 



(18) 



The probability distribution of the number of observed clus- 
ters A^oba is Poissonian, and can be simplified as before in 
our working limit of one or zero observed clusters per cell. 
The posterior probability distribution was analysed us- 



ing the Laplace's Demon package Hall (20121 for the R sta- 



tistical computing environment ( R Development Core Team 
2012). An adaptive Metropolis Markov Chain Monte Carlo 



algorithm was used, and the resulting constraints on the 
model parameters are given in table [I] The Schechter func- 
tion shape parameters are all consistent with the results 
from the REFLEX data alone, indicating that the WARPS 
data do not provide much extra information to constrain 
those parameters. The confidence contours for the evolution 
parameters are plotted in figure [9] As expected, marginalis- 
ing over the uncertainties on the local XLF reduces the pre- 
cision of the evolution measurements, though the presence of 
evolution (i.e. a difference from A,B = (0,0)) is significant 
at more than 95%. This is the first time that evolution in 
the cluster XLF has included this source of uncertainty. The 
best fitting evolving Schechter function is compared with the 
low-redshift REFLEX Schechter function in figure [10] 

The different evolution models are plotted in figure [TTj 
which compares the number of clusters as a function of red- 
shift predicted by the different model XLFs with that ob- 
served. The no-evolution REFLEX model clearly over pre- 
dicts the numbers of clusters at z > 0.6. There is some 




Figure 9. Likelihood contours for the evolution parameters A 
and B based on a Bayesian analysis and marginalising over the 
uncertainty on the shape parameters of the local XLF. Light grey 
contours show the constraints from the maximum likelihood anal- 
ysis. Contours show the la, 2a, and 3cr confidence limits. 
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Figure 10. Schechter functions with best-fitting evolution pa- 
rameters from the Bayesian analysis. The lines show the form of 
the Schechter function at redshifts 0, 0.5 and 1, and the grey line 
shows the REFLEX Schechter function. 



tension between the z > 0.3 maximum likelihood fit and the 
Bayesian model, driven by the Bayesian model's attempts 
to fit the excess of WARPS clusters at 0.1 < z < 0.3. 



4.2 Excess Number Density in Low-z XLF 

Figures [3] [5j and [IT] show that the WARPS detects a sig- 
nificant excess of systems in the range £x(o.5-2.ok c v) = 
1.0 - 2.0 x 10 43 erg s" 1 and 0.1 < z < 0.3 relative to the 
REFLEX Schechter function. These luminosities correspond 
to ~ 2keV systems, so are poor clusters or galaxy groups. 
Interestingly, a significant excess is seen at the same lumi- 
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0.0 0.2 0.4 0.6 0.8 1.0 



I REFLEX, no evolution 

■ ML evolution 

■ Bayesian evolution 




1e+43 1en 
L„(0.5-2.0 keV) tV-,1 10* 4 erg s" 1 



Figure 11. The observed redshift distribution of the WARPS 
clusters is compared to the distribution predicted by different 
models for the XLF. Both the ML and Bayesian fits are to the 
full sample (z > 0.02). The error bars on the observed counts are 
computed according to|Gehrels[ I f 986 1. 



Figure 12. The observed luminosity distribution of the WARPS 
clusters at z < 0.3 is compared to the distribution predicted by 
different models for the XLF. The error bars on the observed 
counts are computed according to|Gehrels| jl986||. 



nosity range in the 160SD low-z XLF. In determining this 
excess, we have been comparing observed number counts 
to the best fitting model to the observed REFLEX XLF. 
It is worth considering whether an excess is present in the 



REFLEX data, but inspection of figure 4 in (Mullis et al 



20041 reveals that this is not the case (note their figure is 
for H = 50 km s" 1 Mpc" 1 ). 

The best-fitting Bayesian XLF model predicts a sig- 
nificantly larger number of clusters at 2 < 0.3 than the 
REFLEX model, as shown in figure [T2"| This somewhat re- 
duces the significance of the excess clusters in the bump at 
£x(o.5-2.okcV) = 1-0 — 2.0 x 10 43 erg s _1 . Marginalised over 
the uncertainties on the model parameters, the Bayesian 
model predicts 16.5 ± 3.7 clusters, compared with the 28 
observed. As the observed number counts are a Poissonian 
realisation of the model prediction, the probability of observ- 
ing ^ A^obs clusters for a set of model parameters 9 which 
predicts a number A" mo d is 



P(N obs \6) 



P{> AU s |AUd)P(AT mod |0)dA mod . (19) 



We model the first probability distribution as a Poisson dis- 
tribution and the second as a Gaussian with mean 16.5 and 
standard deviation 3.7, which results in a probability of ob- 
serving at least 28 clusters in this luminosity bump of 2%. 
There is thus suggestive evidence that a real excess remains, 
and we now consider possible factors that could contribute 
to this. 

Could the excess be simply a result of sample (cosmic) 
variance? The survey areas of WARPS and 160SD are much 
smaller than REFLEX, but the volumes surveyed are still 
significant. The clusters in the bump feature occupy the 
redshift range 0.1 < z < 0.3 and the volume surveyed by 
WARPS in this redshift range is ~ W r h? Mpc 3 . This vol- 
ume is larger than the volumes in which sample variance is 



expected to be significant (e.g Driver & Robotham 20101. 
Furthermore, only 40% of the WARPS fields and 35% of the 
WARPS clusters are in common with 160SD. This further 



reduces the possibility that the excess seen in both surveys 
is due to sample variance. 

The three dimensional distribution of the clusters in the 
luminosity and redshift range of the excess were examined, 
but there was no evidence for clustering in volume, so the 
excess is not caused by a superstructure of clusters. This 
is expected, since the selected fields are scattered across a 
large fraction of the sky. 

A further possibility to explain the excess numbers is 
contamination in the detected flux from unresolved X-ray 
point sources, for example low luminosity Active Galactic 
Nuclei (AGN). Such contamination was removed where pos- 
sible, but not all contaminating sources are resolved in the 
ROSAT PSPC images. Hence, some residual contamination 
is expected, enhancing the estimated cluster luminosities. 
Detailed modelling of the AGN population is beyond the 
scope of this paper, but this contamination would differ from 
the scatter models discussed above, as the effect is purely ad- 
ditive. Some mass or redshift dependence of the AGN con- 
tamination may be required to manifest the localised excess 
of clusters in the Lx , z plane. 

Finally, we consider if selection bias may be responsible 
for the excess clusters seen in this region of the Lx , z plane. 
This is plausible, given that the excess is close to the flux 
limit in the region of the Lx , z plane where the WARPS is 
most sensitive (see figure [TJ. Eddington bias enhances clus- 
ter number counts when scatter is present in the luminosities 
of the population. The slope of the XLF means that for a 
given flux limit, there are more clusters below the flux limit 
that may scatter into the sample than above the flux limit 
that may scatter out of the sample. There are two sources 
of scatter that may be important: statistical scatter due to 
the counting statistics on the Lx measurement, and intrin- 
sic scatter in the cluster population. Our Bayesian analysis 
allows us to investigate each of these sources of scatter in 
turn, by modifying the model for the population scatter in 
equation \17) . 

Note that if the population scatter is constant with Lx 
and z, then a bias is present at all redshifts and luminosi- 
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ties, and increases towards higher luminosities due to the 
steepening of the XLF, so would not produce a localised 
excess in the L x ,z plane. However, the statistical scatter 
decreases with increasing Lx , as ~ errors reduce above 
the flux limit, and increase below it. This gives rise to a bias 
that is strongest near the flux limit, with counting statis- 
tics allowing clusters nominally too faint for the sample to 
appear above the flux limit, so could plausibly contribute 
to the observed excess. However, the typical measurement 
error on luminosities close to the flux limit is « 15%, and 
even modelling the increasing flux scatter below the flux 
limit, the bias due to this source of scatter was found to 
contribute < 1 additional cluster in the regions of the ex- 
cess, compared to models with no scatter. We note however, 
that our model for the scaling of the statistical scatter with 
flux is simplistic - the measurement errors also depend on 
the exposure time and background level in the source field. 
Modelling this for clusters below the flux limit would require 
extensive simulations and is beyond the scope of this work. 

The intrinsic scatter in luminosity of the cluster pop- 
ulation is known to be significant, and if this varies with 
mass or redshift, it could result in a bias that contributed to 
a localised excess in number counts. The variation of cluster 
scatter with mass is not well measured, but there is evi- 
dence that the intrinsic luminosity scatter decreases above 
z w 0.4 (Maughan 20071, albeit measured in a heteroge- 



neous sample. We test the effect that evolving scatter could 
have on the measured XLF by replacing the scatter model 



in equation ( 17 1 with a lognormal distribution, with a stan- 



dard deviation of 50% at z < 0.45, decreasing smoothly to 
20% at z > 0.55 ( |Maughan||2007| . The effect of this evolv- 
ing scatter is to increase the number of z < 0.3 clusters 
predicted by the Bayesian model by » 10%, but this still 
leaves a significant excess of observed clusters at around 
Lx ( o.5-2.okcV) = 1.0 - 2.0 x 10 43 erg eT 1 . 

Similarly, introducing an ad hoc model of intrinsic scat- 
ter which decreases from 90% to 20% at Lx(o.5-2.okeV) = 
2.5 x 10 43 erg s _1 (designed to maximise the Eddington bias 
effect), only serves to increase the model prediction by 2.5 
clusters in the bins contributing to the excess. When com- 
pared to the REFLEX XLF, the observed excess of WARPS 
clusters remains highly significant in the face of all bias con- 
tributions. 



4.3 Sensitivity to High-Redshift Cool Core 
Clusters 

Since the detection of clusters in X-ray surveys is driven 
by their surface brightness, the presence or absence of cen- 
trally peaked emission associated with a cool core may have 
a strong impact on the detection of a given cluster. There is 
some debate in the literature as to whether there are signif- 
icant numbers of cool core clusters in the high-z (z > 0.5) 
Universe. Both Vikhlinin et al. (2007) and Santos et al 



(20101 find a lack of high-z cool core clusters in the dis- 



tant 400SD sample, but Santos et al. (20101 find evidence 



for moderate cool cores in high-z clusters detected in the 
WARPS and RDCS. They argue that the 400SD (and by 
extension 160SD) may have discarded high-z cool core clus- 
ters as being unresolved, but note that this does not imply 
incompleteness in the 400SD providing the surface bright- 
ness dependence was modelled into the selection function. 



The excellent agreement of the evolution in the XLF 
seen in the WARPS and 160SD at z > 0.3 despite their 
very different cluster detection algorithms implies that this 
is indeed the case, and that there are no significant prob- 
lems with the selection function of either survey. We also 
note that the comparisons of clusters detected or missed by 
WARPS and 160SD in 157 common fields revealed no signif- 
icant discrepancies, with differences being explained by the 
differing selection criteria ( jHorner et al.|2008| ). 



5 CONCLUSIONS 

We measured the evolution of the XLF out to z ~ 1 from 
the combined WARPS-I and WARPS-H surveys, finding sig- 
nificant evidence for negative evolution, in the sense of a re- 
duction in the number density of massive luminous clusters 
with redshift. This is confirmed by comparing expected and 
observed numbers, and more convincingly by the maximum 
likelihood analysis. This is consistent with previous measure- 
ments of the evolution of the XLF, and the expectations of 
hierarchical structure formation in a ACDM Universe. 

We investigate the sensitivity of these results to vari- 
ous sources of systematic uncertainty affecting the WARPS 
XLF and selection function. The results are not significantly 
affected by the modelling of the core radii of the clusters, the 
assumed upper redshift limit of the survey, or the technique 
used for estimating detection volumes for the clusters. The 
assumed value of z max — 1.1 is fairly conservative, as there 
is a cluster detected at z — 1.05 and many of the clusters 
would be detectable in X-rays to significantly higher red- 
shift. Thus the true evolution could be somewhat stronger 
than we measure. 

For the first time, we fully incorporate the uncertainties 
on the low-redshift XLF into a Bayesian analysis of the evo- 
lution, and find that while the precision of the measurements 
is reduced, evolution is still significant at the 95% level. 

The good agreement of the measured evolution in the 
WARPS, 160SD and other surveys suggest that the result 
is not sensitive to the details of the cluster detection and 
follow-up strategy, and that the selection functions of both 
surveys are accurately modelled, including the effects of cool 
cores on the detectability of clusters. 

We identified a significant excess of ~ 2keV systems at 
z < 0.3 observed in both the WARPS and 160SD surveys 
relative to the REFLEX XLF. A Bayesian fit to the WARPS 
data, which uses the REFLEX measurements as priors yields 
a model with slightly higher <f>* and a values than REFLEX, 
and reduces the excess. However, even with this model, and 
with including possible contributions from Eddington bias, 
the excess remains significant. The cause of the excess in not 
clear at present, although its presence in both the WARPS 
and 160SD argue against it being a result of a mis-calibrated 
selection function. New, more sensitive measurements of the 
XLF with surveys like XXL ( |Pierre et aL|2011[ ) and XCS 
( jMehrtens et al.|2011 I will provide better statistics for this 
part of the Lx,z plane, providing a means to address this 
question further. 
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